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ABSTRACT 



We present first results on the observation of 
1809 kcV gamma-ray line emission from the Cygnus 
X region with the SPI imaging spectrometer. Our 
analysis is based on data from the performance ver- 
ification phase of the INTEGRAL instruments and 
comprises 1.3 Ms of exposure time. We observe a 
1809 keV line flux of (7.3 ± 0.9) x 10~ 5 ph cnrV 1 
from a region delimited by galactic longitudes 73°- 
93° and |6| < 7° at a significance level of 8a. The 
1809 keV line appears moderately broadened, with 
an intrinsic FWHM of 3.3 ± 1.3 keV. Although this 
broadening is only marginal (at the 2a level our data 
are compatible with an unbroadened line), it could 
reflect the 26 A1 ejecta kinematics. 
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is produced in this region during hydrostatic nu- 
cleosynthesis in O-type and Wolf-Rayet stars; ex- 
plosive nucleosynthesis in core collapse supernovae 
seems to play only a minor role. Yet, the mea- 
sured 1809 keV intensity is by at least a factor of 
2 larger th an expected from curre nt nucleosynthe- 
sis models ijKnodlseder et afll2002|) . This may hint 
at an enhanced efficiency of hydrostatic 26 A1 nucle- 
osynthesis, such as expected for rotati onally induced 
mixing in fast rotating O-type stars l|Vuissoz et all 
2004). Any nucleosynthesis process that is at work 
in Cygnus should in principle also hold for the entire 
Galaxy. Hence Cygnus X may turn out to become 
the Rosetta stone for understanding 26 A1 nucleosyn- 
thesis galaxy wide. It is therefore particularly impor- 
tant to accurately measure the 1809 keV gamma-ray 
line emission from this region, to precisely determine 
its intensity, to identify emission counterparts, and 
to study the line profile that holds information about 
26 A1 ejecta dynamics. 



1. INTRODUCTION 



The Cygnus X region is one of the most active nearby 
star forming regions in our Galaxy (Kn odlsedeit 
12003!) . It houses an exceptionally large number of O- 
type stars that enrich the interstellar medium with 
fresh nucleosynthesis products, either through their 
strong stellar winds, or eventually through their fi- 
nal supernova explosions. The radioactive isotope 
26 A1 is a key tracer of this process, and the study of 
its 1809 keV decay radiation is a powerful probe for 
nucleosynthesis in this region. 

COMPTEL measurements ha ve revealed th e pres- 
ence of 26 A1 i n Cygnus X (pel Rio et all 1 19961 
iPluschkel l2001j) . A deep study of the stellar popu- 
lations in the Cygnus field suggests that most 26 A1 



2. ANALYSED DATA 



The data that we analysed in this work were taken 
during the INTEGRAL performance verification 
phase. They span the INTEGRAL orbits 14-25 (ex- 
cluding orbit 24), covering 1.3 Ms exposure time. We 
analyse single-detector and double-detector data (SE 
and ME2). Energy calibration has been performed 
orbit-wise, resulting in a relative (orbit-to-orbit) cal- 
ibration precision of ~ 0.01 keV and an absolute ac- 
curacy of ~ 0.05 keV (Lonjou et al., these proceed- 
ings). 
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Figure 1. SPI single- detector spectrum in the area around 1809 keV together with the background model, 
(black/red histogram: line component; light/green line: continuum component; grey/red histogram: summed 
model). The shaded areas indicate the energy bands that have been used for modelling the continuum emission 
underlying the line. Instrumental background lines have been labelled by their source isotopes. 



3. BACKGROUND MODELLING 



The instrumental background of the SPI telescope at 
1809 keV is composed of a broad line complex on top 
of a flat continuum distribution (cf. Fig.QJ. The 1809 
keV complex is a blend of at least 3 individual lines: 
a yet unidentified component at about 1805.5 ± 0.5 
keV, a line at the 26 A1 rest energy of 1808.65 keV at- 
tributed to 26 Na(/3~) 26 Mg decay and to excitation of 
26 Mg nuclei, and a line at 1810.77 keV attributed to 
56 Mn(/T) 56 Fe decay l|Weidenspointner et all [2003(1 . 
The relative contributions of each of the lines to the 
complex amount to roughly 10%, 60% and 30%, re- 
spectively, and appear stable with time. 

We model the energy and time dependence of the 
instrumental background using a two component 
model that allows for different time variabilities of 
line and continuum background. The line compo- 
nent is predicted from empty field observations that 
are assumed to be free of celestial 1809 keV emission. 
For this purpose, a total empty field exposure of 4.1 
Ms has been collected by gathering all data for which 
SPI pointed at galactic latitudes \b\ > 20° (covering 
INTEGRAL orbits 14-106). The spectra for this ob- 
servation are summed for all detectors, and the line 
component is extracted by subtracting a constant off- 
set whose level has been estimated from the adjacent 
energy bands 1786-1802 keV and 1815-1828 keV. The 
time dependent counting rate of the line component 
has been modelled by adjusting a multi-component 
background variation template to the counting rate 
history of the empty field data. The components of 
the template reflect our current knowledge about the 
physical processes that give rise to the 1809 keV line 
complex (see also Knodlseder et al., these proceed- 
ings). 

The continuum component was assumed to be con- 
stant in energy (note that we also used a constant 
to extract the line component from the empty field 
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Figure 2. Residual SE and ME2 SPI energy spec- 
tra after subtraction of the background model. For 
the SE spectrum a degradation correction has been 
applied (section^. 



data, hence the line and continuum components are 
complementary and the sum of both should add 
without bias). Its time dependent counting rate has 
been estimated by adjusting the germanium satu- 
rated event rate as activity tracer to the counting 
rate in adjacent energy bands (the same intervals 
have been used as for the extraction of the line com- 
ponent from the empty field data). The adjustment 
has been performed using an adaptive running av- 
erage that requires at least 10000 counts in the ad- 
jacent bands. In this way, short term variabilities 
are explained by the activity tracer while long term 
variations come from the data themselves. 

Figure shows the background subtracted spec- 
tra for single-detector and double-detector events. 
Clearly, a residual signal remains at 1809 keV in the 
spectra of both event types. 
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Figure 3. SE residual spectrum after degradation correction, based on a \ 2 minimisation of the residuals in the 
energy interval 1758-1770 keV (this energy band is indicated by a shaded area). 



4. DEGRADATION CORRECTION 



The exposure of the SPI camera to high-energy par- 
ticle radiation in space leads to a degradation of the 
energy resolution with time due to cre ation of trap- 
ping sites in the germanium crystals ijLeleux et all 
2003). These defects are regularly remov ed using an 
annealing procedure ijRoaues et all 12003). Since the 
empty field observations used for modelling the shape 
of the instrumental 1809 keV line complex were taken 
at difference epochs than the source data, and hence 
at different detector degradations, the shape of the 
instrumental line differs slightly between the source 
data and the background model. A close look to 
Fig. H reveals indeed that the background model 
shows broader lines than the source data (detector 
degradation leads to an extension of the left wing of 
a line due to incomplete charge collection). 

To match the spectral shape of source data and back- 
ground model, we therefore reassign energies E — > E' 
for each event using a random number generator that 
follows the probability density function 



P(E'\E) = 



-(E-E')/AE 



(1) 



where E' < E (i.e. incomplete charge collection) and 
AE characterises the detector degradation. Since 
Eq. ^ effectively decreases the energy of each event 
(by AE on average) the background model has to be 
shifted to lower energies in order to match the data. 
Also, a recalibration of the data is required after this 
correction, that we perform using the 1764.36 keV 
background line of 205 Bi. 

The degradation constant AE and the background 
shift were determined for each detector by minimis- 
ing the residuals in the 1764.36 keV background line. 
This has been done using a \ 2 minimisation ap- 
plied to the energy interval 1758-1770 keV. The re- 
sulting average degradation correction amounted to 
AE « 0.5 keV, the average shift of the background 



model was —0.3 keV. The corrected SE residual spec- 
trum is shown in Fig. Note that the background 
model components were scaled by the fitting algo- 
rithm for x 2 minimisation and visual display since 
the normalisation of the 1764.36 keV line differs from 
those of the other background lines due to a differ- 
ent time history (without scaling the spectrum could 
never become flat). For this reason a strong residual 
remains in Fig. at 1779 keV and also the resid- 
ual 1809 keV line signal is exaggerated. For data 
analysis, however, no scaling has been applied to the 
background model. 

Remaining residuals of the 1764.36 keV line are of 
the order of 1% of the instrumental background line. 
The expected celestial 1809 keV signal amounts to 
~ 12% of the counts in the 1809 keV complex. As- 
suming that the background subtraction leaves also 
1% residuals at 1809 keV, we estimate residuals from 
instrumental background that amount to ~ 10% of 
the celestial 1809 keV signal. 



5. SPECTRAL ANALYSIS 



Straight background spectra subtraction does not 
make use of the coded- mask and field-of-view proper- 
ties of SPI; on the other hand, image deconvolution is 
feasible only for strong signals. As we expect a weak 
signal, we derive spectra by fitting the amplitude of 
a given skymap of (expected) intensity distribution, 
per 1 keV wide energy bin. As skymap we used the 
Dirbe 240 /zm intensity distribution, since this map 
has been shown to provid e a reasonable tracer o f 
galactic 1809 keV emission l|Knddlseder et all ll9991. 

Figure 01 shows the resulting SE and ME2 spectra. 
Fitting Gaussian shaped lines on top of a constant 
provides flux, line position and line width estimates 
that are summarised in Tabled The results for SE 
and ME2 have been combined by weighting with the 
statistical uncertainties. The line position and width 
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Figure 4- SPI SE and ME2 spectra of the Cygnus X region assuming an 26 Al line intensity distribution that 
follows the Dirbe 240 /im infrared emission as a tracer of the massive star population. 



Table 1. Spectral analysis results. The flux is given 
in units of 10~ 5 ph cm ,s~ l and has been deter- 
mined by integrating over galactic longitudes 70° -93° 
and laltitudes b < 7°. The energy and intrinsic line 
width are given in units of ke V. 



Events 


Flux 


Energy 


FWHM 


SE 


6.7 ±2.3 


1808.4 ±0.4 


3.2 ± 1.7 


ME2 


8.1 ±3.1 


1808.4 ±0.7 


3.4 ±2.4 


Sum 


7.2 ±1.8 


1808.4 ±0.3 


3.3 ±1.3 



values are given after degradation correction and re- 
calibration. The line width is the intrinsic value af- 
ter (quadratic) subtraction of the instrumental line 
widths. The instrumental widths have been deter- 
mined from the width of the 1764.36 keV background 
line (after degradation correction and recalibration) 
which is known to be intrinsically narrow. We obtain 
instrumental FWHMs of 3.2 and 3.5 keV for SE and 
ME2, respectively. 



6. DISCUSSION AND CONCLUSIONS 



The measured 1809 keV line flux corresponds well 
to previous COMPTEL measurements. From two 
years of COMPTEL data, iDel Rio et all l)1996ft de- 
rived a flux of (7.0 ± 1.4) x 10~ 5 ph cm _2 s _1 from a 
sky region covering I = [73°, 93°] and b = [— 7°,7°], 
corresponding to a detection significance of 5c Inte- 
grating the fitted Dirbe 240 fim map over the same 
area results in a SPI value of (7.2 ±1.8) x 10~ 5 
ph cm _2 s _1 . Our quoted flux uncertainty also in- 
cludes the uncertainty in the line width. Neglect- 
ing the line width uncertainty by analysing SPI 
data in single energy bands leads to a flux value of 
(7.3 ± 0.9) x 10~ 5 ph cm~ 2 s^ 1 , corresponding to a 
SPI detection significance of 8a. 



From the analy sis of 9 years of COMPTEL data, 
iPluschkei ll200l derived a flux of (10.3 ± 2.0) x lO" 5 
ph cm _2 s _1 from a slightly larger sky region covering 
I = [70°, 96°] and b = [-9°, 25°]. For the same re- 
gion we find a value of (10.0±2.6) x 10~ 5 ph cm _2 s _1 
from the spectral analysis, and (10.3 ±1.3) x 10~ 5 
from the single energy band analysis. Thus, our flux 
measurements are consistent with earlier COMPTEL 
results, and although only 1.3 Ms of exposure time 
have been analysed so far, they present the most pre- 
cise measurements of the 1809 keV line intensity from 
Cygnus X that is available today. 

SPI measured for the first time the precise energy 
of the 1809 keV line towards Cygnus, and the re- 
sulting value of 1808.4 ± 0.3 keV is compatible with 
the rest energy of the 26 A1 line at 1808.65 keV. 
When interpreted kincmatically, the SPI measure- 
ment corresponds to a radial velocity of —41 ± 
50 km s -1 which is compatible with the radial ve- 
locity of ~ km s" 1 of massiv e star clusters found 
in Cygnus X ^Kuodlsederl l20rB. 




FWHM (keV) 

Figure 5. \ 2 of the combined SE and ME2 spectral 
fit as function of intrinsic 1809 keV line width. 

SPI also provides for the first time an estimate of 
the intrinsic width of the 1809 keV line in Cygnus X. 
Figure El shows the x 2 of the combined fit of SE and 
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ME2 spectra as function of the intrinsic line width. 
No significant (i.e. more than 3a) line broadening can 
be claimed from the present data, yet there are hints 
for a modest broadening. The formal line width of 
3.3 ± 1.3 keV corresponds to a Doppler broadening 
of 550 ± 210 km s" 1 (FWHM). This translates into 
expansion velocities of 170-380 km s _1 if a thin ex- 
panding shell is assumed, or to 240-550 km s _1 for a 
homologously expanding bubble. 

These values are relatively high with respect to ex- 
pansion velocities of a few tens of km s _1 that would 
have been expecte d from wind blown bubble expan- 
sion in Cygnus X l|Knodlsedeit 12003!) . Galactic ro- 
tation can also not explain this broadening since the 
velocity dispersion in this area of the Galaxy is quite 
small (also only a few tens of km s" 1 ). Yet turbulent 
motions in hot superbubbles can reach velocities of 
a few 100 km s _1 (De Avillez, private communica- 
tion), and the 26 A1 ejecta may eventually follow these 
motions. However, more observations of Cygnus X 
are required to confirm this hypothesis, by providing 
more stringent informations on the 26 A1 line profile. 
More observations of Cygnus X have been executed 
during the INTEGRAL AO-1 cycle, and further ob- 
servations are scheduled for AO-2. In the near future 
it can therefore be expected that SPI will provide 
a more constraining measure of the 26 A1 line width, 
providing valuable information on the 26 A1 ejecta dy- 
namics and the ISM kinematics in the Cygnus X re- 
gion. 
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